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I. INTRODUCTION 



We apply a random matrix model to the study of the phase diagram of QCD with two colors, two 
flavors, and a small quark mass. Although the effects of temperature are only included schematically, 
this model reproduces most of the ground state predictions of chiral perturbation theory and also 
gives a qualitative picture of the phase diagram at all temperatures. It leads, however, to an 
unphysical behavior of the chiral order parameter and the baryon density in vacuum and does 
not support diquark condensation at arbitrarily high densities. A better treatment of temperature 
dependence leads to correct vacuum and small temperature properties. We compare our results 
, at both high and low densities with the results of microscopic calculations using the Nambu-Jona- 

■ Lasinio model and discuss the effects of large momentum scales on the variations of condensation 
^SJ ' fields with chemical potential. 

D 

A number of early and recent model calculations indicate that quark-quark correlations may play an important 
PsJ role in QCD at finite density. By inducing pairing gaps A ^ 100 MeV, such correlations may prove relevant for the 
J> , physics of neutron and compact stars as well as that of heavy-ion collisions (2|-||] . Due to the limitations of standard 

■ Monte Carlo techniques when applied to systems with finite baryon density, the study of these effects through lattice 
\ simulations is very challenging. A non-zero chemical potential leads to a complex determinant of the Euclidean Dirac 
r operator and results in massive cancellations among configurations. This difficult problem stands in the way of a 

qualitative understanding of the critical physics involved in lattice simulations . 
^-H . One response to this situation is to consider QCD-like theories with additional antiunitary symmetries that guaran- 
' tee the reality of the fermion determinant and therefore render such theories more amenable to lattice formulations. 
\ Examples of such theories include QCD with an arbitrary number of colors, N^, and adjoint quarks as well as QCD 
^ with two colors and fundamental quarks. Knowledge of the critical physics in these cases may offer clues regarding 
I ' the critical physics of the more difficult problem of QCD with three colors. 

Here, we will consider QCD with two colors. In this case, quark and antiquark states transform similarly under 
' global color rotations. They can be combined into spinors with an extended fiavor symmetry SU(2Afy) for which {qq) 
• • ■ baryons and {qq) mesons belong to the same multiplets [ |To| , pd] |. In particular, the lightest baryons and the pions have 
' ^ common mass, to^. This spectrum determines the properties of the ground state for small chemical potential /i > 0. 
, General arguments p2[ indicate that a transition from the vacuum to a state with finite baryon density should take 
5^ place at a critical chemical potential, /ic, which is the lowest ratio of energy to baryon number that can be realized 
^ \ by an excited state of the system. Here, this state is populated by light {qq) baryons, and the transition is thus 
expected to take place at /j,c = TOtt/S. The (T, /x) phase diagram of QCD with two colors has been studied previously 
by Dagotto et al. using a mean-field model of the lattice action [Q. Their results confirm the transition at jic- More 
recently, the smallness of /ic ~ m7r/2 has been exploited in studies of the zero temperature phase transition using 
chiral perturbation theory extended to the flavor symmetry SU(2A^/) |l^JT^. Many of these model calculations have 
been verified by recent lattice simulations [p^ - [2l| . 

We have recently constructed a random matrix model of QCD with three colors which permits thermodynamic 
competition between chiral and diquark condensation p2| . The extension of this model to QCD with two colors 
is straightforward and constitutes a useful consistency check. The random matrix model respects the global flavor 
and color symmetries of the QCD interactions and eliminates much of their remaining detailed structure. For two 
colors, this model offers a qualitatively correct picture of the phase diagram but, at flrst sight, appears to fail in two 
respects. First, in it simplest implementation in which only the two lowest Matsubara frequencies are retained, we 
find an unphysical behavior of the baryon density and the chiral condensate in vacuum. Second, the random matrix 
interactions saturate for large fi and do not support diquark condensation at arbitrarily high baryon density. 

Our primary goal in this study is to examine these difficulties and their relative importance in establishing a reliable 
picture of the phase diagram. The first difficulty (i.e. the ground state behavior of the baryon density and the chiral 
condensate) can be resolved by including all quark Matsubara frequencies, which is equivalent to the introduction 
of proper Fermi occupation factors. For small T, the full frequency sum also modifies the temperature dependence 
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of the onset chemical potential from a quadratic law with fJ,c{T) — Mc(0) ~ to essential singular behavior with 
fJ.c(T) — fj.c(0) exp(— S/r), where E is the vacuum chiral condensate. The remaining phase diagram is qualitatively 
unchanged. The Matsubara sum is equivalent to a monotonic mapping of the phase diagram of the original random 
matrix model, a mapping which therefore preserves the topology of the phase structure. 

The second difficulty (i.e. the saturation of the interactions at large /i) results from the lack of a true Fermi sea 
of states in random matrix models. We consider the effects of a proper Fermi sea by studying a microscopic model 
with an interaction similar to that of the Nambu-Jona-Lasinio (NJL) model. This model naturally exhibits diquark 
condensation at all densities. It also yields an additional power law correction to the essential singular dependence of 
/Ltc(r). It will be shown that other deviations from the prediction of the random matrix model arise only in regions 
where the condensates are weak and thus sensitive to the detailed form of the interaction. Also related to non- universal 
features of the interactions is the observation that deviations from the results of chiral perturbation theory for /i on 
the order of /ic are determined by the form of the interaction on momentum scales substantially larger than fic- This 
is true for both random matrix and NJL models. The sign of these deviations depends on model details and cannot 
in general be established from first principles. 

The remainder of this paper is organized as follows. We review the random matrix model of Refs. extend it to 
the case of two colors, and discuss its low temperature properties in Sec. |l[ We analyze the sum over all Matsubara 
frequencies in Sec. |l| and study the NJL model in Sec. [lV[ We discuss the phase diagrams and the condensation 
fields for each case and present our conclusions in Sec. M. 



II. THE RANDOM MATRIX MODEL 

We first consider the random matrix model of Refs. |2|] in a form appropriate for two colors and two flavors. In this 
model, the single-quark Lagrangian is represented by a matrix whose block structure reflects the global chiral and color 
symmetries of QCD. The matrix elements describe the degrees of freedom associated with the background of gluon 
fields. Aside from the constraints imposed by the prescribed block structure, all matrix elements are independent and 
chosen at random. An exact solution can be obtained if their distribution is Gaussian. Then, an integration over this 
distribution produces a partition function of the form 

Z{fi,T) = J da dA exp [-A^r2(cr, A)] , (1) 

where T is the temperature and fi is the chemical potential. Here, N is proportional to the matrix size and represents 
the number of low-lying degrees of freedom involved in the phase transition. This number scales with the volume of 
the physical system, and the thermodynamic limit therefore corresponds to oo. The integrals in Eq. (|l|) are 

performed over auxiliary chiral and diquark fields, a and A. We consider here the familiar chiral channel, a ~ (99), 
and the diquark channel, A ^ (9 C75 ■''2 o), where C is the charge conjugation matrix {Cjf^C ~ 7J) and where 
T2 and T2 are the antisymmetric SU(2) color and flavor generators ]l5| . (For a related random matrix model which 
takes into account additional effects from instanton-anti-instanton molecules, see p3| .) 
For single-gluon exchange, the thermodynamic potential il assumes the form 

n{a, A) = A{a^ + A^) - Tr log S{a, A) + ^cg, (2) 

where A characterizes the coupling strengths of the interactions in the chiral and diquark channels. (These strengths 
are equal for single-gluon exchange.) The quantity S{cr, A) is the propagator for a single quark in the background of 
the condensation fields, cr and A. The trace is performed over the quark Matsubara frequencies tOn = i^J.+ (2n + l)7rT, 
where n is an integer running from —00 to 00. 

The random matrix model describes only the thermodynamic contribution from low-energy modes. The effects of 
the remaining high-energy modes can, in principle, lead to a "regular" term, fl^eg, which does not affect the critical 
physics. This term is independent of a and A and is an analytic function of T and fi. Although the form of llreg 
cannot be determined by random matrix arguments, we retain it as a reminder of the presence of high-energy modes. 
We will later exploit this degree of freedom to adjust the small temperature properties of the model. 

In its simplest implementation, the model retains only the two lowest Matsubara frequencies n — ±1. The trace in 
Eq. (^ is then reduced to 

Tr log S{a, A) = ^ \og{El {a, A) + n^T^), (3) 
± 

where E±{a, A) represent the quark (— ) and antiquark {+) excitations energies in the background of the condensation 
fields. 
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£;±(cr,A) = {{a + m±fif + A 



2\l/2 



(4) 



The simplest model thus consists of Eqs. (||), (H) and (^). We note that a similar form of the potential fl{<T, A) is also 
obtained for a chiral random matrix theory with a Dyson index P = 1 , in which one retains only the chiral and 
diquark channels. In contrast to the present model, this theory makes explicit use of the pseudoreality of QCD with 
two colors by choosing real matrix elements. However, it does not refer explicitly to color and spin quantum numbers. 
As a result, the diquark condensate is different from that in the present model and develops instead in the ((7^759) 
channel. Since the two models possess interactions with similar symmetries among their respective condensates, they 
naturally lead to similar potentials ^^(ct. A). 

Returning to our model, we first consider the chiral limit, to = 0. Then, the condensation fields in vacuum 
(/i = T = 0) appear only in the combination cr^ + A^, and the system obeys the extended flavor symmetry SU(4). 
The spontaneous formation of a chiral condensate, a and A = 0, further breaks the symmetry down to Sp(4). A 
pure chiral solution, (ct. A) = (S,0), can be transformed into a pure diquark solution, (a. A) — (0,S), by an SU(4) 
rotation. These two solutions are thus thermodynamically indistinguishable, and they describe a single phase with 
cr^ + A^ = = 2 /A. This phase persists up to a critical temperature Tc = S/tt above which all condensation fields 
vanish and the SU(4) symmetry is restored via a second-order transition. (This is a mean- field result; renormalization 
group arguments show that for two or more flavors, fluctuations may actually drive the transition first-order p5| .) 
For T — 0, sl positive chemical potential, > 0, breaks the SU(4) symmetry down to SUl(2) x SUr(2) x U(1)b; a 
quark mass m ^ breaks it further to SUv(2) x U(1)b [^^^- In the following, we consider light quarks by choosing 
TO <C E, where E = ■\/2/A is the vacuum chiral condensate for m = 0. Note that for m 7^ 0, the chiral field is no 
longer a true order parameter, and the second-order phase transition at /i = becomes a cross-over. 

The diquark field. A, remains a true order parameter even when to 7^ 0. The system favors the formation of diquark 
pairs in regions where d^fl/dA^{a,0) < 0, and a pure chiral solution can be a saddle-point of the potential at best. 
There is therefore a second-order phase transition along the line d^il/ dA^ {(t(p,T) , 0) = 0, where a{fi,T) obeys the 
chiral gap equation dil{a, 0)/da = 0. This line is given by the relation 
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and delimits a region in the {T,fi) plane as illustrated in Fig. 1 (a). (We will discuss Figs. 1 (b) and (c) below.) Its 
intercept with the T axis determines the onset chemical potential 
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which becomes fJ.c{T) ~ /ic(l -I- 7r'^T^/(4E^)) for small T. According to the arguments in the introduction, the onset 
chemical potential of Eq. (^) also obeys Hc = to^/2. This leads us to identify the pion mass as to^ ~ to^/^E^/^-|-C'(to). 

This result is not trivial and may be regarded as the random matrix analogue of the current-algebra relation 
™^ ~ TO^/^. To understand the origin of this relation in the present model, consider the random matrix vacuum. For 
simplicity, assume one quark flavor and ignore diquark condensation by setting A = 0. Consider first to = 0. In this 
case, we know that the spontaneous breaking of chiral symmetry guarantees the existence of a massless Goldstone 
mode [ p6| independent of the number of colors | p7[ |. Consider next the case of a finite but small mass to, for which 
chiral symmetry is explicitly broken. The previous Goldstone mode now describes the lowest-lying fluctuations around 
the global minimum of Sl(cr, 0), which is given by (Jo = S - to/2 + 0{m^/E). From Eqs. (|), (|), and (|) with A = 0, 
we find that the energy of fluctuations around this solution is given as 



fluct 



2 

E2 



(((To + Serf + Sn'^) - 2 log (cto + m + Serf + Sn 



(7) 



where the pion mode St: describes oscillations in the direction with (qj^q) 7^ and 5a describes those in the (qq) 
direction. Thus the mass, to, acts as an external field with an almost linear coupling to the order parameter; for small 
TO, the symmetry breaking term goes as (Jfifluct cx; — toctq- This coupling causes the curvature of the free-energy in 
the soft direction to be proportional to to. Expanding to second order in da and Sn and working to lowest order in 
m/T,, we find 



n\ 



Sa—S-TT—O 



E2 



Sa' 



2m 
E3 



(8) 



This allows us to identify the mass ratio m^/m^ as ~ to/E and deduce toJ ^ to, which confirms the result of Eq. (|^). 
The generalization of this argument to include fluctuations in the {qq) channel is straightforward and leads to results 
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similar to those of chiral perturbation theory for the constant terms ||T^. (Random matrix theory cannot reproduce 
terms with spatial derivatives.) 

We now return to the diquark condensed phase. The two gap equations, dfl/da — and dfl/dA — 0, possess a 
single solution with both 7^ and A 7^ 0, 

cr = -TO+— — (9 

2 /i^ — 

where I] = (2/j4)^/2 Expanding to the lowest orders in m/E, the fields near the condensation edge /i « /Zc obey 

CT^ + A^ = s2-7r2T2 + 0(mE). (12) 

For r = and to leading order in m/S, these results coincide with chiral perturbation theory ||l5| . For /i > fic, the 
chiral field decreases like and this saddle-point solution rotates into a pure diquark solution keeping cr^ + A^ 

nearly constant. Corrections of order mS in Eq. ( p^ depend on the explicit functional form of the potential of 
Eq. (y) (i.e. here, a logarithm involving certain combinations of a and A) and are thus not universal. We will 
see below examples of other functional dependences which give different corrections in mE. If we now turn to the 
non-condensed phase (/i < /ic and A = 0), we find 

TO M2_^2y2 fm^\ 

where terms of order m and higher are again not universal. 

Although the results to lowest order in m/T, are consistent with our expectations, higher-order terms give rise to 
certain unphysical properties. For T = and /i < /ic, we expect that, in the absence of confinement, the lowest-lying 
excitations carrying a net baryon number are single quarks with an energy S. Below the condensation edge, we have 
fj, < ^ and the excitation energy ~ E is well above the Fermi level /i. The ground state should thus have a chiral 
field equal to its value at /i = and a baryon density of zero. We find instead a chiral field cr ~ E — to/2 -|- /i^/ (2S), 
which increases quadratically with fj,. The behavior of the baryon density, ne = —dQ/dfi, is also incorrect. From the 
solutions above, we find that the potential in each phase obeys 

4to ^2 - TT^T^ 

L 2 — 

to2 u'^+Tt'^T'^ /to2\ 
n ~ f^vac --^-^ ^^^2 + ^reg + O i—j {fl > fie), (15) 

where the vacuum energy is fJvac = 2 — 21ogE2. These relations lead to negative baryon density for /i < /ic! 

This behavior is a result of having truncated the Matsubara sum in Eq. (^ to the two lowest frequencies. This 
approximation does not affect the critical physics but docs modify bulk properties. It is therefore appropriate to 
correct the shortcomings of such an approximation naturally with an appropriate adjustment of the regular term, 
ilrcg, in the thermodynamic potential of Eq. (^). The requirement that the baryon density should vanish for /i < /tc 
is readily met by the choice flrcg = —2[i}lY?. For T = 0, we then obtain a baryon density of 

n^-J-^-[ + O(ip^) (/.</.c) ,^g. 

which has the same /t-dependence as that found in chiral perturbation theory |l5[] . For /i < //c, the terms of order 
©(TO/ic/E^) are positive, hence riB > for all [i. 

To summarize, the random matrix model reproduces the results of chiral perturbation theory to leading order in 
m/E for the condensation fields and for the baryon density, provided only that we exploit the freedom of including a 
component f^rcg in the thermodynamic potential. This addition cannot, however, affect the condensation fields, and 
a still grows weakly with /t below the condensation edge. The variation of the condensation fields near the edge is 
shown in Fig. 2 (a). (Figs. 2 (6) and (c) will be discussed below.) It is, in fact, difficult to discern the dependence of 
(7 on /i for /t < /tc as the quadratic term ~ A<^/S is of order to/E relative to the leading term, which is of order ~ E. 
Figure 3 (a) shows the condensation fields for selected temperatures. 



17 ~ f^vac - — + 2 + (/i</,c), (14) 
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III. SUMMING OVER ALL MATSUBARA FREQUENCIES 



The /^-dependence of a and found in the preceeding section was obtained by retaining only the lowest Matsubara 
frequencies. Although this approximation is reasonable for the study of the finite temperature phase diagram of QCD 
with three colors, it is bound to break down in the present case for small temperatures near the condensation edge. 
As T becomes smaller, adjacent Matsubara frequencies ±37rT, ±57rT, . . . come closer to those retained and eventually 
form a continuum as T ^ 0. In physical terms, the frequency continuum at T = enforces the exclusion principle by 
introducing Fermi occupation factors. 

As shown in the Appendix, the sum over all Matsubara frequencies leads to a partition function of the form 

Z{^i,T) = J dadA exp{-Npn{a,A)} (17) 
where (3 = l/T. The thermodynamic potential is now given by 

r!(a. A) = ^^^^ - {^±('^' ^) + '^T\og[l + exp{-pE±ia, A)}]} + fi.eg • (18) 

The excitation energies, E±{a, A), are given by Eq. (Q). As before, we allow for the inclusion of a term fl^-cg, analytic 
in /i and T, which represents the contribution from high-energy modes. The requirement that the T = /i = vacuum 
has zero baryon density is met by taking Slj-cg as a constant. This constant is further set to f2rcg = — 2m, in order 
for the vacuum to have zero pressure. Again, E is the chiral field for T = ^ = m — 0. 

Compared to the previous random matrix model, the full Matsubara sum has introduced Fermi occupation factors, 
f{x) = {1 -|- exp(2;/r)} ^ . This is easily seen in the two gap equations 

In this form, the gap equations possess solutions which, at T = and for all n below the condensation edge, exhibit 
the same properties as those of the vacuum, /i = T = 0. For /i < /Zc, we now find a pure chiral solution with 
A = and a constant chiral field cr = S. The baryon density vanishes exactly. The onset chemical potential, 
/ic, at T — is determined by requiring that both gap equations are satisfied while taking A — > 0^. We find 
— im{I] + m)y/^ ~ (mE)^/^ -|- 0{m). For fj, > the gap equations possess a single solution which reproduces 
the result of chiral perturbation theory to lowest order in m/E. Higher-order terms depend on the functional form 
of Eq. (^8|) and are different from those obtained in Eqs. (^) and (|l^). These terms are therefore not universal. We 
show the fields a and A near the condensation edge in Fig. 2(6). Figure 3(b) presents their behavior as a function of 
fi for selected temperatures. 

The phase diagram in the full (T^fi) plane, shown in Fig. 1(b), appears as a mere remapping of the diagram of 
Fig. 1(a). Both diagrams have the same topology. The most dramatic effect of this remapping is a sharp vertical 
stretching of the condensation edges at low and high chemical potentials. We show in the Appendix that, at the lower 
edge, the previous behavior of fic{T) — Atc(O) ~ is now replaced by essential singular behavior with Hc{T) — /Xc(0) ^ 
exp[— (I]-|-m)/T] sinh[/ic(0)/r]. This behavior is clearly a consequence of the introduction of Fermi occupation factors. 
These factors automatically imply detailed balance, whereby two single-quark states separated by an energy difference 
AE are weighted by a relative factor of exp(— A£'/T). Except for these modifications at low temperature, the global 
picture remains remarkably similar to that of Fig. 1(a). 



IV. BUILDING A FERMI SEA 



Cooper pairing occurs in any attractive channel independent of its interaction strength. Thus, even though the 
QCD interactions weaken on the momentum scales appropriate for large densities, the existence of a diquark phase 
at mean field level should naturally continue to arbitrarily high For three colors, model calculations with large /j, 
and small coupling constant, g, find a pairing gap A ~ fig~^ eyip{—c/g) where c is a constant [ p8[ . Similar behavior 
is to be anticipated in QCD with two colors. 
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By contrast, the random matrix interactions saturate at large and do not support diquark gaps at arbitrarily 
high density. This weakening of (qq) correlations is ultimately related to the absence of a Fermi surface. The quark 
states in the random matrix model do not carry a momentum, and there is therefore no Fermi sea. In fact, for fixed 
a and A, the system has only two energy levels (the E±{a, A) in Eq. (^), whereas a microscopic model exhibits a 
continuum of states labeled by their particle momenta. 

In this section, we study the effects of a true Fermi sea on the global picture established above and consider 
a microscopic theory with an interaction suggested by the Nambu-Jona-Lasinio model p9| ]. The structure of the 
interaction for large momenta must be specified by an appropriate form factor. A sharp three-momentum cutoff 
saturates the interactions at large momenta and does not lead to a diquark gap for all fi. Since we wish to explore a 
model which supports pairing at all densities, we choose instead a soft form factor, 

which weakens the interactions with a scale of g ^ A. This choice leads to a model similar to that adopted by Berges 
and Rajagopal in their study of diquark condensation for ^ 3 Our intent is to investigate general trends, 

and we are not concerned with precise parameter values. Thus, we take A = 3S, where S is again the chiral field in 
vacuum for m = 0. The mass, m, is chosen as m = 0.01 E in order to have a clear separation between the three scales 
TO, VtoS, and S. (We have explored other ranges of parameters and found the same basic trends.) The introduction 
of momentum dependence transforms the model of Eq. ( p^ into 

n ^ A{a^ + A^) - I dqq^ {E±{q, a, A)+2T\og{l + exp{^f3E±{q, a, A)}) j, (22) 

where the excitations energies are now momentum dependent, 

E±{q, a, A) - {{E{q) ± m)2 + T\q)A^) , (23) 

with E{q) — {q^ + (to + crJF2((7))^)^/^. As before, the coupling constant A is fixed by requiring that the vacuum chiral 
field be equal to cr = S for to = 0. 

The system now naturally supports a diquark condensate at all densities. To illustrate the mechanism at work, 
consider first the chiral limit to = 0. Then, for T = and /i > 0, we find that cr — while A obeys the gap equation 

A = dqq^YTjqf , . (24) 

When A = 0, the right side of this equation diverges logarithmically due to the behavior of the integrand for q near 
/J,; it vanishes for large A. Hence for A > 0, the equation always has a solution with A > 0. For small /i, A is 
sufficiently large that the logarithmic singularity for q ^ ii is inoperative. The range of contributing momenta is 
then determined by the form factor J-^{q) and includes all g < A. As increases, the momentum window in which 
the integrand is logarithmically large, |g — ~ AJ^^(/i), decreases due to the form factor, and the strength of the 
logarithmic singularity increases. It eventually dominates the integral at large /i and leads to exponentially small gaps 
which decrease like A(^) cx exp[— ^/(^^J^(/i)^)]. 

A similar mechanism is encountered when we depart from the chiral limit. For to > 0, we have a ^ 0, and 
the condensation fields are given by the solution of two coupled gap equations. These equations again exhibit a 
logarithmic divergence as A ^ for momenta q satisfying E{q) ~ where E{q) — {q^ + (to + crF(g)^)^)-'^/^. We have 
^'('Z) ?^ ™ + where the chiral field is of order a ~ S for jj. ^ jic and decreases for larger ^. Hence, the condition 
E{q) ^ /i is not obeyed for small fi as E{q) E when — > 0. It is, however, satisfied starting at some /i in the 
range fic < fJ' < ^ where a decreases below E. For large fi, the singularity condition can always be met. Therefore, 
for T = and /i > E, the logarithmic singularity is active and dominates the momentum integral. It follows that 
A is always non- vanishing for moderate and large /x and that the diquark phase persists for all large fi as shown in 
Fig. 1(c). 

Figure 3(c) shows the condensation fields for various fixed T. As we have just argued, the maximum diquark 
field, and therefore also the critical temperature, decrease exponentially with fi when /i > A. In this region, the 
thermodynamic potential is dominated by large momenta and is very sensitive to form factors. This situtation stands 
in contrast to the three color case , which has different coupling constants and exhibits a different evolution of the 
diquark field. There, most of the interesting variations in A(/i) and in the phase transition line take place below the 
cutoff, /i < A. These results are likely to be much less sensitive to form factor effects. 
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Returning to the problem at hand, it is useful to remember that the present models have been implemented at 
mean-field level. Results in regions where condensation fields are weak cannot be regarded as reliable since weak 
condensates may not survive fluctuations. An illustrative example of the effects of long-wavelength fluctuations has 
recently been considered in lattice studies of the NJL model in two dimensions. There, the data suggests a strong 
suppression of the diquark order parameter with respect to the mean-field predictions and also provides indications 
of strong critical fluctuations |^ . 

Results for small /i are more reliable. The critical physics near the condensation edge, fj,c, is primarily determined 
by symmetries and is most likely protected by them. For the NJL model, we find results which are remarkably similar 
to those of the previous sections. For zero T, we find a second-order phase transition from a chiral phase to a diquark 
phase at a critical chemical potential /Zc ~ rn^l'^ . The chiral field is constant for T = and fi < fic and decreases 
roughly like 1/^^ for fi > fi^. The diquark field and the baryon density increase in this region. For T = and near 
the condensation edge, the fields follow chiral perturbation theory as illustrated in Fig. 2 (c). For small T, the critical 
chemical potential behaves as ^c(r) — ^c(O) ^ T^/^ exp {— S/T} sinh[/Xc(0)/T]. The exponent 3/2 in the power law 
correction is directly related to the dimension of the momentum integrals in Eq. (^) . 

V. CONCLUDING REMARKS 

The two random matrix models and the NJL model considered here all reproduce most the predictions of chiral 
perturbation theory in the vicinity of T = and near the condensation edge, fi — fic- However, the matrix models 
are essentially different from chiral perturbation theory on two respects. First, matrix models are mean field and 
do not treat the dynamics of long wavelength Goldstone excitations. Unlike chiral perturbation theory, the matrix 
models are thus not sensitive to effects related to spatial derivatives of the Goldstone fields. Second, matrix models 
are formulated on a large range of fi and T which extends far beyond the condensation edge. As a result, they contain 
additional, non-universal corrections to chiral perturbation theory. These corrections can, in part, be attributed to 
the effects of momentum scales larger than m^. 

This can be seen by considering both chiral limit and departures from it. In the strict chiral limit, m = 0, the 
condensation edge goes to /ic = 0, and much of the critical structure which we have previously discussed disappears. 
The chiral field vanishes for all /i > 0, and a finite diquark field develops for T below a critical temperature Tc{fJ-). 
Variations in A(/i) as a function of fi occur on a large and model-dependent scale. Away from the chiral limit, we 
can distinguish two regions. First, in the vicinity of the condensation edge, all thermodynamic quantities follow the 
predictions of chiral perturbation theory and vary with fi on a scale fic ~ TOtt- Second, far from the condensation edge, 
the thermodynamic properties of all three models vary on the much larger scales encountered in the chiral limit. In 
fact, for fi ^ fic the condensation fields differ from those for m = by powers of m. The scale ~ m^/^, which is 
central to chiral perturbation theory, thus no longer plays a role in this region where larger and non-universal scales 
are dominant. 

The resulting picture for fi ^ fic found for each of the models considered here is thus quite different from that 
predicted by chiral perturbation theory. For the random matrix models, the relevant scale of variation is the vacuum 
chiral field S, which is intimately related to the variance of the matrix elements. This is also the scale on which (qq) 
correlations weaken. The saturation of the interactions manifests itself in a decrease of A(/i) with fi in both random 
matrix models, compare Figs. 2 (a) and (b). By contrast, the dominant scale of variation for /x 3> /Xc in the NJL 
model is the momentum cutoff, A. This scales acts in a different way. Because of the logarithmic singularity in the 
gap equations discussed above, A(/i) increases until /i ~ A. For larger values of /i, interactions weaken due to the 
form factors, and A decreases with fi. The logarithmic singularity also causes the baryon density to rise as fi increases 
from fic- For the NJL model, deviations from the results of chiral perturbation theory are opposite to those of the 
random matrix models as shown in Fig. 1(c). 

The lesson is that large momentum scales tend to play a dominant role away from the condensation edge. This 
should be kept in mind when interpreting lattice data. In this regard, it is interesting to note that the lattice mean- 
field action of Dagotto et al. leads to a phase diagram with a second-order line which terminates at some high 
fi and T = in a manner similar to the phase diagrams in Figs. 1(b) and (c). This ending of the transition line is 
attributed in Ref. ||l^ to a saturation of the lattice interactions. 

A final comment concerns the evolution of the baryon density as a function of fi. We have not displayed baryon 
densities for fi ^ fic because random matrix models describe only the contribution of soft modes. Since the baryon 
density receives contributions from both high- and low-energy modes, its behavior can be quite different from the 
predictions of random matrix models. Similarly, baryon densities calculated with the NJL model are sensitive both to 
the form factor chosen and to variations in the quark mass m over much of the phase diagram. Again, this suggests 
that such results should be considered with care. By contrast, we have demonstrated that thermodynamic properties 
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near the condensation edge are dominated by symmetries and that the models considered provide a quahtatively 
correct description of the baryon density in this region. 

In conclusion, we have found that a random matrix model which retains only two Matsubara frequencies correctly 
reproduces the leading-order results of chiral perturbation theory. Beyond leading order, however, it produces certain 
unphysical results including negative baryon densities for small fi and a chiral condensate which varies with fi. These 
pathologies can be eliminated by the inclusion of all Matsubara frequencies. This extension of the random matrix 
model also correctly produces an essential singularity in the second-order line at T = 0. Neither variant of the 
random matrix model considered here can support diquark condensation at arbitrarily high baryon density since 
neither contains a true Fermi sea of states. A microscopic model capable of supporting diquark condensation at 
all densities does so at the cost of introducing significant model dependence. Near the condensation edge, a three 
dimensional model also induces power law corrections to the temperature dependence of the onset chemical potential. 
This result is robust. 

The present random matrix models contain the ingredients necessary for a mean-field description of the critical 
physics near the condensation edge. All thermodynamic quantities in this region are dominated by those symmetries 
which we have implemented in the interactions. Given only a minor and physically reasonable adjustment in the 
treatment of the Matsubara sum, the model reproduces the results of chiral perturbation theory or of any other 
microscopic model which implements the same symmetries. Away from the condensation edge, however, the random 
matrix interactions saturate and lead to non-physical results. It seems unlikely that "improvements" of the theory 
in this region, whether in the form of an extended random matrix model or of a detailed microscopic model, will be 
beneficial. The large n region is not dominated by symmetries, and results in that region will unavoidably be fragile 
and model dependent. Since the random matrix approach is by construction free of the inevitable details of more 
microscopic models, it provides a simple and useful way of distinguishing between those results which are dominantly 
influenced by symmetries and those which are not. 



A recent study of SU(2) lattice gauge theory with four quark flavors by J. B. Kogut, D. Toublan, and D. K. 
Sinclair [ p3| has revealed the existence of a tricritical point at which the phase transition changes from second to 
first order as temperature is increased. An approach based on the Landau free energy for the diquark condensed 
phase shows that such point is allowed by the symmetries. It arises from a term which contains odd powers of the 
order parameter and which is related to the Goldstone fields that carry the same quantum numbers as the diquark 
condensate. 

The mean field level of the present random matrix models does not include the dynamics of the Goldstone fields. 
Thus, our models do not naturally lead to a tricritical point. When interpreted as a Landau free energy, the expression 
for the thermodynamical potential, r^(i7. A), contains only even powers of A and hence precludes the presence of such 
a tricritical point in the phase diagram. 



We thank D. Toublan, J. J. M. Verbaarschot, and A. Schafer for stimulating discussions and K. Splittorff for 
his critical reading of the manuscript. We also thank the DOE Institute for Nuclear Theory at the University of 
Washington for its hospitality and the Department of Energy for partial support during the early stages of this work. 



Here we consider certain technical details of the random matrix model in which all Matsubara frequencies are 
included. Following Ref. |^^, we introduce a Euclidean time which is the Fourier conjugate of the Matsubara frequency 
and which runs from < t < P = 1/T. 

Consider first the logarithmic term in il{a,A), Eq. (0), whose frequency trace is 
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APPENDIX 



A. Sum over all Matsubara frequencies 



oo 




log[/32 ((2n + i fiT^T^ + A2 + (fj + m ± fif)] . 



(25) 



zb n— — oo 
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Here, we have included a piece of the term fiiog from Eq. (||) by inserting the prefactor = in order to ensure 
that the argument of the logarithm is dimensionless. The sum over n on the right side of this equation is best evaluated 
by taking the derivative of TrlogS' with respect to either cr or A. Using then the summation formula 



oo 



t--hN= ,2 + (, + 1/2)2,2 ^ (26) 

n— — oo 

and integrating back over the condensation fields, we arrive at the result 

Trlog S{a, A) = {£;±(a, A) + 2riog[l + exp[-/3S±(a, A)]]} + c(m, T). (27) 
± 

Here, c is independent of a and A and the E±{a, A) are the single-quark energies of Eq. (^). 

The frequency sum over the quadratic terms is simpler. Since the condensation fields are constant in time, Fourier 
transformation to Euclidean time gives 

TiAia^ + A^)^ dtA(a^ + A^) ^ pA{a^ + A^). (28) 







The saddle-point effective potential is the sum of Eqs. ( p7| ) and (|2^) plus an additional regular piece, /3f2rog, which 
includes c(/i, T) and is chosen as follows. We wish to obtain a T — fi — vacuum to have zero baryon density. This 
constraint is met by taking fij-eg as a constant. The additional requirement that the T = fi = vacuum should have 
zero pressure sets firog = — S — 2m, where E is the vacuum chiral field. Since we are primarily concerned with the 
functional form of the effective thermodynamic potential, we have chosen to tune A to A = 1/S in order to obtain a 
chiral field cr = E in vacuum. This choice simplifies comparison with the other models. Removing a common prefactor 
(3, we arrive at Eqs. ( p^ ) and (|l8|). 



B. Temperature dependence of the onset chemical potential 

In this section, we determine the temperature dependence of the onset chemical potential, ^c{T), for both the 
random matrix model with all Matsubara frequencies included and the NJL model. We first illustrate the derivation 
for the simple random matrix model. For fixed T, the onset chemical potential satisfies the chiral and the diquark 
gap equations with A set to zero. From il(cr. A) in Eq. (|l8|), we have 

- a = tanh[ + tanh[ , (29) 

- ^{a + m + ^l)-Han\l[ + (a + m - /i)"' tanh[ 1^^]- (30) 

The small T behavior can be obtained by expanding the chiral field and the chemical potential around their values 
at T = and ^ = jic (i.e. cr = S + (5(T and = + 5^) and expanding the explicit temperature dependence as 
tanh[a;/(2T)] ~ 1 — 2exp[— x/T]. To first order, this expansion gives 

2^ ~ -4exp[-(cr + TO)/T]cosh[^c/r], (31) 



~ -5a 



1 1 



(E + TO + ^c)^ (E + m — ^c)^ 
1 1 



(E + m + /ic)^ (E + TO — /ic) 
2 exp[-(E + TO)/T] 1 1*:^ . (32) 

Inserting the first of these relations into the second and expanding the denominators to first order in /ic ^ m}/'^ and 
zeroth order in to, we obtain the result given in the text: 

<5A.~Eexp[-^]sinh[^]. (33) 
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Due to the sensitivity of the NJL model to both the mass m and the cutoff parameter A, it is more difficult to 
obtain an exact relationship for ndT) in this case. We can, however, get a rough estimate by following the previous 
steps and restricting our attention to terms of zeroth order in m/E and first order in /ic/S. An expansion of the two 
gap equations to these orders then gives 

5^i ^ csinh[^] dqq^^ eM^E{q)/T] , (34) 

where c is the constant c = {J^ dqq'^!F{q)^/E{q)^)~^ and E{q) = (g^ + E^J^(g)^)^/^ is the single-quark energy for 
A = and m = 0. For T — > 0, the integral on the right side can safely be approximated by saddle-point methods. 
Expanding E{q) as E + q^(l — 4E^/A^)/(2E) and setting q = in the denominator, we obtain 

. Ai, (T\"^ E exp[-E/r] 

Comparison with Eq. (|3^) reveals an additional power-law correction of T^/"^ . This correction is a consequence of 
having a genuine continuum of momentum states; the exponent 3/2 is half the dimensionality of the space. 
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FIG. 1. Phase diagrams in the (T, /i) plane for the random matrix model with two Matsubara frequencies included (a), with 
all frequencies included (b), and for the NJL model (c). In each case, the diquark-condensed phase with A 7^ is separated 
from that with A = by a second-order line. In Figs. 1 (a) and (6) , E is the vacuum chiral field in the limit m = 0. In Fig. 1 
(c), A is the momentum cutoff used in the NJL model. In each of the three cases, To is the critical temperature for /u = and 
in the chiral limit m = 0. 
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FIG. 2. Condensation fields near the edge (T, ^) = (0, /ic) in the random matrix models with two Matsubara frequencies (a), 
all Matsubara frequencies (6), and in the NJL model (c). In each case, the dashed lines represent the predictions from chiral 
perturbation theory and the continuous lines are the results from the models considered in this paper. In each panel, ao is the 
chiral field at /u = 0. The grey curves represent the densities as a function of /tt, scaled so that their slope is unity at the edge 
fj, = fic- Note that the density curves from chiral perturbation theory overlap with those from the random matrix models ((a) 
and (6)). 
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(a) 




fi/A n/A 

(e) (/) 

FIG. 3. Chiral and diquark fields for selected temperatures in the random matrix model with two Matsubara frequencies 
included ((o) and (6)), with all Matsubara frequencies included ((c) and (d)), and in the NJL model ((e) and (/)). In a given 
graph, the outer curve corresponds to T = 0. Going inward, the inner curves correspond to T/To — 0.2,0.4,0.6,0.8, 1.0, 1.2, 
where To is the critical temperature for ^ = and m = 0. Note that the a curves slightly overlap in (c) for T/To — and 
T/To = 0.2 and those in (e) slightly overlap for T/To = 0,0.2, and T/To = 0.4. In Figs. 3 (6), (d), and (/), the diquark field 
vanishes for a few of the highest selected temperatures. In each graph, S is the vacuum chiral field in the limit m = 0; in 
Figs. 3 (e) and (/), A is the momentum cutoff used in the NJL model. 
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